Magnetic fluctuations and resonant peak in cuprates: a microscopic theory 
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The theory for the dynamical spin susceptibility within the t-J model is developed, as relevant for the resonant 
magnetic peak and normal-state magnetic response in superconducting (SC) cuprates. The analysis is based on 
the equations of motion for spins and the memory-function presentation of magnetic response where the main 
damping of the low-energy spin collective mode comes from the decay into fermionic degrees of freedom. It is 
shown that the damping function at low doping is closely related to the c-axis optical conductivity. The analysis 
reproduces doping-dependent features of the resonant magnetic scattering. 
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Since its discovery in inelastic neutron scattering experi- 
ments in superconducting (SC) YBa 2 Cu 3 07 [|lp, the magnetic 
resonance peak has been the subject of numerous experimen- 
tal investigations as well as theoretical analyses and inter- 
pretations. The magnetic peak has been systematically fol- 
lowed in YBa2Cu306+ K (YBCO) into the underdoped regime 
[||, ||, £|], where the resonant frequency oj r decreases while the 
peak intensity is increasing. Its pronounced appearance is still 
related to the onset of SC, although it could start appearing 
even at T > T c . More recent results confirm similar behavior 
in Bi2212 and T12201 cuprates [§]. 

Several theoretical hypotheses have been considered for 
the origin of the resonant peak: that it is a bound state in 
the electron-hole excitation spectrum [H], a consequence of a 
novel symmetry between antiferromagnetism (AFM) and SC 
|^] and that it represents collective spin-wave-like mode in- 
duced by strong AFM correlations [[| ^]. There is also an on- 
going debate whether the resonant peak is intimately related to 
the mechanism of SC and whether it can account for anoma- 
lies in single-electron properties, as tested in angle resolved 
photoemission spectroscopy. 

The scenario of a resonant mode as a collective magnetic 
mode seems to correspond well to experimental facts, in par- 
ticular the qualitative development of the resonant mode with 
doping and its onset for T < T c . Still the status of the the- 
ory of the resonant mode, and moreover of the magnetic re- 
sponse in cuprates in general, is not satisfactory, both from 
the point of understanding and of the appropriate analytical 
method. Relevant microscopic models, such as the Hub- 
bard model and the t-J model have been so far studied in 
the weak coupling or random-phase approximation |^], ne- 
glecting strong correlations. The latter have been consid- 
ered using a Hubbard-operator technique JTc|], and more re- 
cently within the self-consistent slave-boson approach JTT|], 
self-consistent spin-fluctuation method [|l2[, as well as within 
the phenomenological spin-fermion model [|^, 0], 

Our aim is to develop a theory of the dynamical spin sus- 
ceptibility Xq(w) within the t-J model. The natural approach 
to analyse collective modes is the memory-function formalism 
[13]. In analogy to the previous study of spectral functions 
[14] we employ the method of equations of motion (EQM) 



to generate the spin dynamics and in particular to establish 
the effective decay of localized spins into fermionic degrees 
of freedom, which is the essential ingredient to describe the 
damping of the collective mode as well as the destruction of 
the long-range AFM order at finite doping. 
We study the t-J model 
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on a square lattice. In order to be close to the physics in 
cuprates we take into account besides the nearest neighbor 
(n.n.) hopping tij — t also the next n.n. fy = t' hopping. 
Strong correlations among electrons are incorporated via pro- 
jected operators, e.g. c\ s = (1 — a )ct a , which do not 
allow for double occupancy of sites, and in J <C t where for 
cuprates we assume furtheron J = 0.3t. 

Within the memory function approach of Mori [ pj| to dy- 
namical response functions, which we here pursue, the dy- 
namical spin susceptibility x q (w) can be expressed as 
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in which form it is particularly well suitable for the analysis 
of coherent collective magnetic response, as manifest in the 
resonant peak in cuprates. Here, u q is related to the disper- 
sion of the collective mode provided that the mode damping 
is small enough, i.e. 7 q ~ M q '(o;q) < w q . In the opposite 
case, we are dealing with an overdamped mode, as seems to 
be generally the case for the magnetic response near the AFM 
wavevector q ~ Q = (tt, tt) in the normal state of cuprates. 

In order to evaluate the quantities entering Eq.(Q) we follow 
the formalism of memory functions Jl3|], defining the scalar 
products and projections in terms of static response functions 
(A\B) = x°ab = -<(^ f ; -S)) w=0 , and the action of the Liou- 
ville super-operator LA = [H, A]. Within this framework we 
can expres 

r /q = (CS^\C SQ = CS$]) , u\ = ry q / Xq , 

Af q H= (Q£ 2 S%\[£q ~ u^lQ^S^) /?7 q , (3) 
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where Cq — QCQ is the projected Liouville super-operator, 
Xq = Xq( w = 0) is the static susceptibility and the projectors 

Q = Q'Q are given by Q = 1 - I^M" 1 ^! and Q' = 
l-l^))^]" 1 ^!, respectively. 

In order to proceed we write down equations of motion 
(EQM) for the spin operators S^. By evaluating CS^ it is 
straightforward to explicitly express ry q in Eq.(|^) as 

V* = Tn Ek> q+ k + £ q _ k - 26 k ](45 ks > 

+ ^E k [^+k + J q -k-2J k ](5+Sl k ), (4) 

where e k is the 'free' band dispersion following from tight- 
binding hopping tij. We are interested in the vicinity q ~ Q. 
We note that ?7 q is closely related to the internal energy, i.e. 
VQ ~ ~~ (H) ■ This indicates that i]q is not strongly dependent 
neither on temperature T nor on q, while the doping depen- 
dence Ch = 1 — (rii) (at low doping) rj ~ ac^jt| + 6 J is as 
well modest, where a~6 ~1, as inferred e.g. from numerical 
studies of the model. 

The evaluation of £ 2 S q is also straightforward, but more 
tedious. It is convenient to separate the action L 2 = C 2 + 
£f + £j with C 2 = [C t ,Cj]+, involving different powers 
of kinetic and exchange terms, respectively. In the following 
analysis we assume that for the damping function r q (u>) = 
M''(uS) at q ~ Q, at low T ~ and in the doping regime of 
interest, i.e. of low to intermediate doping, the essential term 
is C 2 . Namely, we presume that at least in the normal state 
the damping r q (w) approaches a constant 7 q for ui — ► as 
in a Fermi liquid (although anomalous). Such a damping can 
arise only from the coupling to fermionic degrees of freedom, 
which then has to involve H t term. 

One can give more arguments in support of our assumption. 
In an undoped system - AFM, only C 2 is effective. However, 
it is well known that at T = such a term gives a damp- 
ing 7 q oc q 2 where q = q Q, as well as Tq(u}) oc w 2 [[l5|] . 
Hence it leads to underdamped and well defined AFM magnon 
excitations for q — ► 0, the origin being in the phase space of 
low-lying spin excitations restricted to q ~ 0. The argument 
about the C 2 j term can be extended to a doped system. As- 
suming the Fermi liquid form Xq( w ) oc a mode-coupling 
treatment here would again yield r q (cj) oc uj 2 . The role of 
C 2 , on the other hand, cannot be apriori neglected at lu — ► 0. 
However, within the same decoupling described below we get 
a vanishing contribution of diagonal Cf S% to M q (w). 

We should stress that the evaluation of C 2 S* requires ex- 
plicit consideration of the projections of fermionic operators 
in H t . In the site-representation we can express 

A 2 Sj = -E fc ^(s J -- s k )v jk + 

+ J2ki s tjk[tkiSjTjk^j s c ls - tjiTjiB^CkaSk] + H.c. (5) 

Here 7y = m(l — rij) + V%j and Vij = riirij/2 + 2Sj-Sj- 
is the spin-interchange operator. Complicated form of Eq.fph 
reflects the well-known involved nature of correlated hopping 
in a strongly correlated system, i.e. with a reshuffling of spins 
along the hole path. Since the main goal is to get the coupling 



to nonlocal fermionic degrees, in EQM we replace operators 
Vij, %j by their thermal averages Pij, Tij, respectively, lead- 
ing to an effective hopping renormalization in Eq.(||). Note 
that for the n.n. hopping, Ti represents essential reduction, 
since in a Neel state one would get Pi = and only in heav- 
ily doped system Pi ~ (1 — c/j) 2 /2. For the regime of interest, 
i.e., 0.1 < Ch < 0.25 one can on the basis of numerical re- 
sults for the t-J model simplify T\ ~ Ch- We have to apply 
to Eq. (Q) also the projector Q. To the lowest order this im- 
plies that the projected operator does not contain explicitly the 
initial operator S* itself. So we get 

^kq=(ek-4 +q )(et-?k +q )-Cq ; 

where e k is defined with renormalized hopping parameters 
w UjTij whereas £ q is determined by the condition 
Ek Wk q = 0- We note that within an analogous approxi- 
mation CjS* = 0. 

Eq.(§) represents a decay of spin variables into fermions 
in a doped system, so performing a decoupling (in a normal 
state) we get for the damping in the lowest approximation 

r q M = i~ I — [/(" - a;') - f(co')]R^, J) , (7) 
with 

«q(w, W) = ^ ^k(a;')Ak +q (a; - uf) , (8) 

k 

where A k (w) are electron spectral functions. At low doping 
an alternative decoupling of fermionic operators directly in 
the site representation, Eq.(Q), neglecting the coherence be- 
tween different sites might be more appropriate, considering 
the fact that in underdoped systems spectral functions exhibit 
pronounced incoherent behavior. This would yield 

Jtq(w, w') Rj ir(lN{u;'W(u - J) , (9) 

where Af(u>) = (2/iV) Ek-Ak^) is the electron density of 
states (DOS). The form (g, fyis particularly appealing since 
by Eq.(^) the damping becomes proportional to the c-axis con- 
ductivity, i.e. r q (cj) oc u c (uj), which can be well represented 
within the same form [flq]. 

Before proceeding to discussion of r q (cj) we note that the 
theory requires an additional input, i.e. in Eq.(^|) we need 
either u; q or x q - A possibility is to fix unknowns with the sum 
rule in the paramagnetic phase 

i r°° 

~ J o d0J Cth — Xq(^) = (S'_&) = Cq , (10) 

where in addition (l/AOEq^q = 0- ~ c h)/4. Static cor- 
relation functions C q are rather well known within the t-J 
model [pp7|], in particular one can express Cq oc £ 2 where £ is 
the AFM correlation length, determined also experimentally 

in La 2 - K Sr 2; C04 with (, oc 1/^/qT. 
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In the present study we assume some simple forms for 
spectral functions Ak(ui) although there exists also an anal- 
ogous approach to fermion dynamics, yielding as well pseu- 
dogap (PG) features in A^(uj) from the coupling to spin de- 
grees Jl4|]. The simplest assumption in the normal state is 
to insert some effective coherent band crossing the Fermi en- 
ergy, i.e. Ak(oj) ~ Z^6(uj — e k ). Such a form yields at 
T -» r Q (w -» 0) ~ const from R, Eq.(g), as well as 
from expression ([)]), provided that the Fermi surface crosses 
the AFM zone boundary. The doping dependence enters T in 
several ways. First, there should be an overall proportional- 
ity to Ch which is evident from the relation of T with a c . In 
an effective band picture the latter can arise from vanishing 
QP weight Z tx ^Jch assuming that the effective band width 
W ~ ach\t\ + bJ does not vanish for Ch — > 0. r\ and are 
less Ch-dependent in the regime of interest. We can estimate 
the size of normal-state damping as 7q ~ 7t£qZ 2 /2W 2 i]q, 
which at low doping can be quite small 7q <C t. Neverthe- 
less, from available numerical data we estimate 7q still too 
large for an underdamped collective mode at ojq to exist. 

In order to get un underdamped resonant mode at ui r one 
needs a depleted damping TQ(uj r ). The latter can evidently 
arise in the SC state, T < T c , from the SC gap. We can 
model this by introducing into Eqs.(|8|^|) an effective d-wave 
gap Ak = Ao(cosfca; — cosfc y )/2 via Ak(a>) and J\f{u>). 
Due to the broken symmetry we include in Eq. ^ also the 
anomalous spectral functions Fk(u>) [S]. The SC gap even- 
tually leads to the vanishing of Tq(ui < uiq) — where 
Wq ~ 2Ak* < 2Ao while k* is the position of the 'hot spot' 
along the AFM zone boundary. 

Such an analysis with a SC d-wave gap is particularly ap- 
propriate for the situation close to optimum doping. Results 
for r q (w) for several q along the zone diagonal are presented 
in Fig. la. Parameters are chosen so as to correspond to op- 
timum doping Ch ~ 0.2, i.e. effective band with t = 0.3i, 
i = -0.lt, Z = 0.4, and we take A ~ O.lt. We note that 
the damping with a single step at q = Q develops two steps 
for q 7^ Q and the threshold oj* — > closes for q = |q — Q| > 
q* ~ 0.3 [|J. This is a mechanism for the onset of strong 
damping of the collective mode for q > q* and its disappear- 
ance as the resonant character emerges. Another important 
feature is the 'normal state' damping 7q = Tq(uj > oJq) 
calculated in this approach explicitly. We note that 7q 3> uj r , 
preventing any coherent feature in the normal state. 

The corresponding Xq( w ) is presented in Fig. lb. We 
choose w q ~ 0.38t such as to yield C q ~ 0.4, being con- 
sistent with results within the t-J model for intermediate dop- 
ing [|l7|]. At q = Q the resonant mode is undamped since 
2Ao > oj r - It should be stressed, however, that due to quite 
large Tq(ui > luq) the resonant frequency is significantly 
renormalized 



uj r = luq[1 + M'q(u)/u\ 
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while its intensity is also reduced Iq ~ CQUj r /ujQ. The rest 
of the spectral weight is distributed over a shallow but very 
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Figure 1: (a) Damping function r q (aj) at optimal doping c; t ~ 0.2 
for various q || Q for a d-wave SC, (b) corresponding spin response 
Xq(^). Note that t ~ 400meV. The resonance peaks are artificially 
broadened with 8 — O.Olt. 



broad continuum. On the other hand, moving away from Q 
the mode gets overdamped and merges with a broad contin- 
uum for q > q*. In Fig. lb we as well observe a downward 
dispersion of the resonant peak consistent with experiments 

Analysing the regime of low doping it appears more ap- 
propriate to use the incoherent approximation, Eq.(^). It is 
crucial that the normal-state damping 7q also decreases with 
doping, scaling approximately as oc c/j. On the other hand, it 
is rather clear that for underdoped cuprates the experimental 
data in the SC phase cannot be explained with a single gap 
only. Neutron scattering results for Xq( w ) i n me underdoped 
YBCO [||] indicate on the appearance of the resonance at ui r 
at T < T c , possibly even at T > T c [|||. However, in contrast 
to optimum doping the resonant mode at oj r is quite damped. 
At the same time, the pronounced shoulder below the reso- 
nance, i.e. lu c < uj r , appears [Q]. The drop in Xq(w < u> c ) 
can be again interpreted with a coherent SC gap in Tq (uj), 
but with substantially diminished gap 2A C < u> r . Since the 
'normal' damping is still too large, i.e. 7q > u> r , we need 
to assume also the appearance of a pseudogap in DOS below 
u> ~ Up for T < T* with T* > T c . This is well consistent 
with the behavior of a c (uj) in underdoped cuprates pT| ] where 
the pseudogap appears at T < T* . In Fig. 2a we present 
characteristic Fq(cj), calculated at T ~ for Ch ~ 0.1 
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Figure 2: (a) rcj(aj) for underdoped c/, ~ 0.1. Dashed curve: nor- 
mal state regime, T ~ T* ~ 0.05i. Full curve: SC regime. Here 
damping is assumed proportional to a c (uj) with an SC gap added 
below lo c ~ O.Olt. (b) Corresponding Xq( w )- 



with corresponding A ~ O.Oli and calculated ijq ~ 0.25i 
whereas tAf(ui > uj p ) ~ 0.15, as inferred from numerical re- 
sults for the t-J model [p2[|. An additional PG reduction for 
u; < ~ O.li is assumed, consistent with experimental data 
pT| ] together with an estimate of the DOS within the pseu- 
dogap tAf(u> < ojp) ~ 0.1 p2]]. The resulting Xq(w), corre- 
sponding to Cq ~ 1.0, is shown in Fig. 2b. In the normal state 
Xq(w) is overdamped. Still, a substantial part of the sum rule 
is exhausted in the window co < wq. It should be also pointed 
out that in the normal state the role of T > is essential 



via the sum rule (10) leading to u>q(T) shifting with T. From 
Fig. 2b it is evident that also in the SC state several features are 
different when compared to larger doping in Fig. lb: a) The 
(spin) gap shoulder appears for u < lo c below which there 
are no spin excitations, b) The resonant peak is damped even 
for T < T c , but still underdamped. c) The spin response and 
the sum rule for Xq( w ) are nearly exhausted within uo < u> p . 
Since in the underdoped regime 'normal' 7q is also reduced, 
peak position given by Eq.([fT]) is not significantly renormal- 
ized and lo t ~ ojq . Consequently, for the underdamped mode 
using the sum rule, we obtain cj r ~ tj/Cq oc 1 /£ 2 oc Ch- 

To summarize, our analysis of the dynamical spin response 
and resonance peak in cuprates within the memory function 
approach can qualitatively, and at low doping even quantita- 
tively, reproduce the spectra as measured in neutron scattering 



experiments [|TJ ||, ||, EP], The central point of the present the- 
ory is the evaluation of the damping function r q (w), using 
the EQM and considering the decay of spin fluctuation into 
electron-hole excitations as the dominant process. In the nor- 
mal state we obtain a large damping 7q > ojq increasing with 
doping, leading generally to an overdamped AFM collective 
mode. Still it is important to realize that the renormalization 
due to Tij <C 1 in Eq.(EJ) is essential. Namely, without this 
reduction the damping would be much too large and in par- 
ticular it would prevent the matching of the sum rule Eq.dlT 
with any pronounced short range AFM order Cq ^> 1/4. 

Addressing the resonance peak in the SC state, we find that 
in optimally doped samples it can arise only inside the fre- 
quency gap ~ 2Ao, being strongly renormalized in compari- 
son to the characteristic frequency ujq and also reduced in the 
intensity. On the other hand the incoherent part of spin fluctu- 
ations extends over a broad frequency range ~ W and for the 
chosen parameters accounts for ~ 80% of the integrated in- 
tensity J Xq(^)^. In the underdoped (weakly doped) case, 
however, the sum rule at T = is almost completely ex- 
hausted within the peak width. As a consequence we obtain 
in this regime to r oc Ch, quite consistent with experiments 
Moreover, in contrast to optimum doping with a single SC 
gap, in underdoped samples two scales seem to play the role 
in the SC state. The proportionality of damping to a c would 
then indicate a gradual closing of the pseudogap leading to a 
crossover into an overdamped situation. More detailed analy- 
sis will be presented elsewhere. 
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